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Abstract 

The rate of secondary charm-quark-pair production has been measured in 4.4 million hadronic 
decays collected by OPAL. By selecting events with three jets and tagging charmed hadrons in the 
gluon jet candidate using leptons and D*^ mesons, the average number of secondary charm-quark 
pairs per hadronic event is found to be (3.20ib0.21ib0.38)xl0~^. 
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1 Introduction 



The production of secondary heavy quarks from a virtual gluon is commonly referred to as gluon 

splitting. This process is considerably suppressed because both the gluon and the quark jet from 

which it originates must be sufficiently virtual to produce the heavy-quark pair. Nonetheless, these 

events make a significant contribution to heavy quark pair production in e+e- annihilation: e+e~— >qqg, 

g— >QQ, where Q is a bottom or charm quark. These events will be referred to here as g— >cc or g— >bb 

events. This paper describes a measurement of the rate of g^cc at LEP at center-of-mass energies in 

the region of the TP peak. 

The probability of producing a heavy-quark pair from a gluon, per hadronic TP decay, is defined 

as 

^ N(Z"^qqg,g^QQ) 

Sqq n(ZO ^ hadrons) ' ^ ' 

This probability has been calculated in the framework of perturbative QCD to leading order in Og, 
with the resummation of large leading and next-to-leading logarithmic terms to all orders . The 
probabilities for the secondary production of a charm-quark or bottom-quark pair are predicted to be 
in the range (1.35-2.01) xlO~^ for gcc, and (1.75-2.90) x 10^^ for g^. Precise measurements of these 
quantities allow an important comparison with QCD calculations, and also reduce the uncertainty in 
the experimental determination of electroweak variables, such as Rc, the fraction of Z"— >cc events in 
hadronic TP decays. 

The first measurement of gcc was made by OPAL [^, where the contribution from the process 
g— >cc to the inclusive D*^ meson momentum spectrum was obtained by subtracting the contribution 
from TP^cc and Z*^— >bb events where the D*^ was produced from a primary quark. This was followed 
by a second OPAL analysis [P, in which gluon jets were selected, and charmed hadrons in these jets 
were identified by a lepton tag. Combining the two measurements gave gcc= (2.38ib0.48)xl0^^ |^], a 
value which is compatible with the upper range of the theoretical predictions. 

The ALEPH and DELPHI collaborations have recently measured the rate of secondary production 
of b-quarkpairstobegj^^=(2.77±0.42±0.57)xl0-3 [§ and g^^=(2.1±l.l±0.9)xl0-3 g] respectively. 
Both measurements are consistent with theoretical predictions. 

The existing OPAL measurements ||5|, ^ only analysed part of the total TP data sample. In this 
analysis the full sample was used. The precision of this measurement was further improved with a 
more refined data calibration, optimisation of the analysis algorithms, and a more reliable Monte Carlo 
simulation of the OPAL detector. To establish a signature for secondary charm-quark-pair production, 
events containing three jets were selected, and the gluon jet identified. The charm content of the gluon 
jet candidate was analysed by identifying electrons, muons, and D*^ mesons, to give three independent 
measurements of g— >cc. The analysis is presented as follows: Section |2| describes the hadronic sample 
and event simulation; Section |^ discusses the selection of events likely to contain a hard gluon and the 
methods used to select the gluon jet; this is followed by a description of the lepton and D*^ channels 
separately in Sections Q and |5|; the systematic uncertainties for all the tagging schemes are listed in 
Section ^; the paper concludes with a summary in Section |^. 



2 Hadronic Event Selection and Simulation 

We used data collected at LEP by the OPAL detector between 1990 and 1995 in the vicinity of the 
Z*^ peak. Hadronic Z^ decays were selected using the number of charged tracks and the visible energy 
in each event as in Reference [10|. This selection yielded 4.41 million events. The primary vertex of 
the event was reconstructed using the charged tracks in the event and the knowledge of the position 
and spread of the e"'"e~ collision point. 

Monte Carlo events were used to determine the selection efficiency and background levels. The 
selection efficiency was measured in dedicated samples of events containing the g— >cc process. For 
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the lepton analysis, at least one of the charmed hadrons was required to decay semileptonically. The 
D*^ analysis used events where at least one of the secondary charm quarks hadronised to a D*^ 
which decayed via D*"*"— ^D'^vr"'", followed by D"— >K~7r"'"0. For background studies, 9 million 5-flavour 
hadronic decays plus an additional 3.5 million Z^^bb events and 2.5 million Z*^— >cc events were 
generated. All these samples were produced with the JETSET 7.4 Monte Carlo program |11|. The 
heavy quark fragmentation was parametrized by the fragmentation function of Peterson et al. [IT 
and the measured values of the partial widths of into qq were used |13] 
different charmed hadrons at ^/s = lOGeV and ^/s = 91 GeV are consistent [14| 
charmed hadrons produced in Z*^— >cc, g^cc, and in b hadron decays were taken from Reference 
as were the semileptonic branching ratios of charm hadrons. All samples were processed with the 
OPAL detector simulation package |15]. 



The production rates of 
so the mixture of 



3 Jet Selection 

3.1 Jet Reconstruction 

Measurement of secondary charm pair production could in principle be reduced to a charm hadron 
counting experiment in events where the primary quarks were light flavoured, i.e. up, down, or strange 
quarks. However, such a sample is hard to obtain as Z*^— >cc and Z'^^bb events can not be identified 
with 100 % efficiency. Instead, by grouping the particles of an event into jets, charmed hadrons which 
are produced from the primary quarks can more easily be identified and counted as background. 

To identify the process of gluon splitting into a charm-quark pair, we required the event to contain 
exactly three jets. Using simulated events, we investigated the effect of various jet finding algorithms to 



select the three-jet topology. Of the Durham, Geneva, Cone, and JADE-EO algorithms |16| that were 
tested, the JADE-EO recombination scheme with a 7/cut value of 0.05 gave the most significant g^cc 
signal. As can be seen from Figure ||(a), the two secondary charm quarks tend to be contained within 
a single jet, so that in three-jet events, g— >cc events were identified more efficiently than background 
events. 

For the selected three-jet events, the jet energies were calculated, neglecting mass effects, using 
the relation 

sm ipjk + sm ipij + sm tpik 

where Ecm is the center-of-mass energy, -i/^jj is the angle between jets i and j and Ei is the calculated 
energy of jet i. This equation only holds for coplanar events and therefore, if the sum of the angles 
between the jets was smaller than 358°, the event was rejected. These criteria retain 30 % of the data 
sample. From the simulation, the efficiency for g— >cc events to pass the selection is 56 %. 

In the previous OPAL analysis @, only 3.5 million hadronic decays were used, with a t/cut value 
of 0.03, compared to 4.4 million events and a ycut value of 0.05 used here. 



3.2 Gluon Jet Selection 

The gluon jet candidate in the event was then selected. Several selection algorithms were tried, of 
which two were retained, based on their efficiency for selecting the correct jet and for background 
rejection. The first method, already used in Reference |Q, assumes that the lowest energy (LE) jet 
is the gluon jet. The second method, used here for the first time, takes the jet that is most readily 
subdivided into two as the gluon jet candidate (JS), and was found to correctly identify the gluon jet in 
g— >cc events more often. This method is motivated by the fact that the gluon jets of interest contain 
two charmed hadrons, and should therefore show more evidence of jet substructure than primary 
charm jets and most bottom jets which contain only one charmed hadron. The JADE-EO jet finding 

^Charge conjugation is assumed throughout this paper. 
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algorithm was applied to each jet individually, and the gluon jet candidate is taken to be the one with 
the highest value of ycut at which the jet splits into two subjets. Using the new JS method, the correct 
gluon jet is identified in 60 % of three-jet events from the signal process g— >cc, while the corresponding 
number for the LE method is 51 %. 

The performance of the two gluon selection schemes was investigated for both the lepton and the 
D*^ analyses. The LE method performed better than the JS method for the lepton analysis, after 
all the additional background rejection criteria had been applied (section 4.2). The overall efficiency 
of the lepton analysis using the LE method was 66 % of the efficiency using the JS method, but 
with a factor of three less background from primary heavy-quark events, yielding a 15% advantage 
in statistical significance. In contrast, for the D*^ analysis, the JS method outperformed the LE 
method; the efficiency for the D*^ analysis doubled when using the JS method while the heavy quark 
background only increased by a quarter, yielding an improvement in statistical significance of 80%. 
As a results of these studies, the LE method was used for the lepton analysis while the JS method 
was used for the D*^ analysis. 

Figure |l](b) shows the ycut value of the selected jet at the point where that jet splits into two 
subjets, for all three-jet events selected from data and the Monte Carlo simulation. Figure ^(c) shows 
the distribution of the jet energy for quark and gluon jets from the process g^cc, and Figure 0(d) 
shows the ycut value at the jet splitting point for these jets. 



4 The Lepton Analysis 

Using the sample of three-jet events, and taking the gluon jet candidate to be the lowest energy jet 
(LE method), events containing the gluon splitting process, g— >cc, were searched for using a lepton 
tag, which assumes that one of the charm quarks decayed semi-leptonically. 



4.1 Lepton Identification 



Electrons were identified using an artificial neural network |17] while muons were identified using in- 
formation from the muon chambers in association with the tracking chambers, as in The lepton 
identification was limited in polar angle 9, defined as the angle between the lepton candidate and the 
beam axis direction of outward-going electrons. The polar angle of identified electrons was limited 
to |cos0| < 0.8, while identified muons were required to be within |cos0| < 0.9. In addition, the 
momentum of the lepton candidate was required to satisfy 2 < p < 9GeV/c, and have a transverse 
momentum, with respect to the jet axis containing the candidate, below 2.75GeV/c. Electrons con- 
sistent with being produced from photon conversions or from Dalitz decays were rejected with an 



artificial neural network using geometrical and kinematical properties |17]. 

The efficiency for leptons from the g— >cc signal events to pass the lepton identification and the 
momentum selection criteria was 16% for electrons, and 17% for muons. This differs from the pre- 
vious OPAL analysis Q, with the most significant difference being the larger momentum range used. 
Previously muons were only accepted in the narrower range of 3 < p < 6GeV/c. The shift to the 
larger momentum range, together with the improvements in electron identification and the photon 
conversion finder, have resulted in a larger efficiency for identifying leptons from g^cc. 

At this point, the data sample contained 7180 tagged electron candidates and 14631 tagged muon 
candidates, corresponding to a gcc purity of 10 % and 6 % respectively, while the efficiency to detect 
a lepton from the g^cc process that passes all the criteria described above (Sections 2-4.1) was 
approximately 7.5%. The large difference between the number of selected electrons and muons is 
mostly due to contamination of hadrons passing the muon selection criteria, as described in section 
4.3. 
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4.2 Suppression of Jet Mis-assignment Background 

With the selection of lepton tagged events described above, background contributions from events 
with jet mis-assignment were unavoidable: Jet mis-assignment events are defined as those where the 
lowest energy jet contains a lepton from the decay of a primary heavy quark rather than from a heavy 
quark in a gluon jet. The jet mis-assignment background is dominated by Z*^— >bb events, and was 
suppressed by means of a b-tagging algorithm based on reconstructed displaced secondary vertices. 
A neural network with inputs based on decay length significance, vertex multiplicity and invariant 
mass information was used to select vertices with a high probability of coming from b hadron 
decays. Events were rejected if any of the three jets were tagged by the neural network. This procedure 
resulted in a reduction of the jet mis-assignment background by 43 %, which included a 54 % reduction 
in bb events, while retaining 78 % of the g^cc sample. After this cut the data sample contained 5 049 
electron candidates and 11031 muon candidates. 

To further reduce the jet mis-assignment background from non-gluon jets, we explored several 
properties of the lowest energy jet, with respect to the other two. Naively, as a jet containing the 
gluon splitting g^cc contains two charm quarks, one would expect the jet mass and multiplicity to 
be larger for that jet than for the other two jets. Comparison of Monte Carlo samples containing 
the signal, with samples of Z^— >bb and Z"— >cc events in which a lepton coming from the decay of a 
primary b or c quark was assigned by the jet finder to the lowest energy jet showed that background 
rejection through mass and multiplicity cuts was indeed feasible, and the following requirements were 
found to give the best background rejection: 

• Max(Mi, M2)/M3 < 2, where Mj is the mass of the ith. jet, and the jets are ordered by energy, 
with jet 1 having the highest energy. 

• (Ni -|-N2)/N3 < 2.5, where Nj is the track and electromagnetic cluster multiplicity of the ith 
jet, where an electromagnetic cluster was counted only if no charged track was associated with 
it. 

These two selection criteria retained 55 % of the g^cc events, while rejecting 71 % of the Z"— >cc 
background and 67 % of the Z"— >bb background. On application of the jet mis-assignment background 
suppression, in conjunction with the lepton tagging selection criteria described previously, and accept- 
ing only one lepton tag candidate in each event (with priority given to electrons due to the higher 
muon fake rate), the data sample was reduced to 2434 electron candidate events and 4362 muon 
candidate events. At this stage, no further background suppression was done. Rather, the remaining 
backgrounds in the data sample were evaluated using Monte Carlo estimates (see Section 4.3). 

4.3 Estimates of Background Rates 

The determination of the rate of each background source is discussed below. A summary of the 
estimated data sample composition is given in Table 1. 

Jet Mis-assignment Background 

The rate of the jet mis-assignment background was estimated from the 14.5 million hadronic decays 
of the Monte Carlo sample mentioned in Section 2, after applying all the selection cuts of Section 4.2. 
The number of leptons from semileptonic decays of charm and bottom hadrons was determined, 
excluding leptons from g— >cc or g— >bb decays. By scaling to the number of hadronic events in the 
data sample, we estimated this background to contain 1 027ib32 candidates from Z"— >bb events and 
702ib26 candidates from TP^cc events, where these uncertainties are statistical only. 
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Quantity 


Electron channel 


Muon channel 


Observed events 


2 434 


4 362 


Jet mis- assignment (cc) 


342±19 


360±18 


Jet mis- assignment (bb) 


494±22 


533±23 


Residual photon conversion 


630±25 




Lepton mis-identification 


81±9 


2 580±51 


Dalitz decays 


149±12 




g ^ bb 


50±7 


53±7 


Estimated signal 


688±42 


836±59 



Table 1: Summary of selected sample sizes and estimated composition. 



Photon Conversions 

From the Monte Carlo simulation, the photon conversion finder fails to tag (14.3ib0.4) % of the con- 
versions, where the uncertainty is from the systematic uncertainties in the conversion finding effi- 



ciency |17|. These untagged conversion electrons then form a background to the electron tagged 
events. Knowing the efficiency of the conversion finder, the background from untagged conversions 
was estimated from the number of tagged conversions to be 630ib25 events. 

Lepton Mis-identification and Decays in Flight 

To estimate the background from hadrons which were erroneously identified as lepton candidates, 
we used the Monte Carlo sample to determine the probability that a charged track with a given 
momentum, p, and transverse momentum with respect to the direction of the associated jet, pt, 
should be incorrectly identified as a lepton. The number of background leptons in the data was then 
derived by multiplying the number of tracks in the data that passed the selection criteria, excluding 
the lepton identification, by these fake probabilities. In practice the fake probabilities per track are 
estimated in bins of p and pt, and corrected for differences between the Monte Carlo simulation and 



the data, as in |18|. The correction for the difference between the MC and the data introduces a large 
uncertainty on the number of hadrons mis-identified as leptons as described in section 6.1 

The total number of hadrons mis-identified as muons was estimated at 2 580ib51. Decays in flight 
of light hadrons into muons are included in this estimate. The number of hadrons mis-identified as 
electrons was estimated by a similar method to be 81ib9. 

Dalitz Decays of 7r° and t] 

The number of background events from the decay of vr'^ and r] into e"'"e~7 was estimated from Monte 
Carlo simulation and corrected to the known vr*^ and rj multiplicities in Z° decays These contri- 
butions were estimated at 149ibl2 electron candidate events. 

Gluon Splitting: g bb 

The number of events from the process g— >bb that survive the selection criteria was calculated from 
Sbb ' ^had ' c'^, where e*^ is the efficiency for at least one lepton from the process g^bb to survive 
the selection criteria, N^ad is the number of hadronic events, and g^^; — (2-69 it 0.67) x 10~^ is the 
averaged measured value of taken from From Monte Carlo simulation, was found to be 
(1.0 lb 0.1) %, leading to an estimated background of 103ibl0. 
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4.4 Comparison of the Data and Monte Carlo 

Since a major part of the background contribution has been estimated from the Monte Carlo, it was 
crucial that the consistency between the Monte Carlo and the OPAL data be verified. Of particular 
importance was the accuracy in the simulation and estimation of the jet mis-assignment background, 
as the efficiency to tag a heavy quark jet as a g— >cc and the fraction of heavy quark jets found in the 
lowest energy jet were critical to this analysis. 

To check the validity of the jet mis-assignment background composition, we examined the lepton 
yield in three-jet events where the lepton did not originate from the lowest energy jet. Specifically, 
we searched for events with leptons in either of the two highest energy jets which passed all but the 
lepton requirement in the third jet of the selection criteria. The leptons were identified using the same 



criteria as described in Section |4.1| , and the jet-based mass and multiplicity cuts were applied in the 
same fashion as in Section However, the search was restricted to prompt leptons from b hadron 
decay by requiring lepton momenta above 5 GeV/c and transverse momenta above 1.5 GeV/c. The 
resulting yield of such leptons per three-jet event was found to be (1.49 ± 0.02) x lO"^ ^.^e Monte 
Carlo and (1.47 it 0.02) x 10~^ in the data. For the Monte Carlo sample the simulation showed that 
the b-hadron purity of these events was 85 %. 

We also looked at events tagged as Z'^^bb with the b-tagging algorithm described in Section in 
order to compare the jet mis-assignment rate in Z"— >bb events. This was done by searching for a lepton 
in the lowest energy jet of events that have vertices compatible with b-hadron decay. After applying 
all other selection criteria we observed 646 such events in the data. The Monte Carlo prediction scaled 
to the data sample size was found to be 632 events. The Monte Carlo simulation estimated that the 
b hadron purity of these events was approximately 93%. From the agreement of the Monte Carlo 
prediction with the data sample estimates for these two tests, the data/Monte Carlo consistency was 
found to be adequate, and the jet mis-assignment well modelled in the Monte Carlo. 

In addition, comparisons of the Monte Carlo simulation prediction and data distributions for the 
variables critical to the lepton tagging analysis were made, with the results of the comparison shown 
in Figure |2|. Specifically, Figure ||(a) shows the ratio of the maximum mass of the first two jets with 
respect to the third jet. Figure ^(b) shows the ratio of the multiplicity of tracks and unassociated 
clusters in the first two jets with respect to the third jet. Figure §(c) shows the lepton candidate 
momentum spectrum for events passing all the selection criteria, and ^(d) shows the lepton candidate 
transverse momentum. All plots present the Monte Carlo prediction which include the g— >cc contri- 
bution normalised to the number of signal events obtained in this analysis. This comparison shows 
an agreement between Monte Carlo and data in both shape and rate prediction. Thus, the procedure 
of subtracting the Monte Carlo prediction for the jet mis-assignment background was justified. 

4.5 Results 

The charm-quark-pair production rate per hadronic event is related to the measured quantities by 

Nscl 



Nhad • e • 2 • B(c ^ Xeu) ' 



(3) 



where the following notation is used; Ngei is the number of events passing the selection criteria after 
subtraction of background events (Table 1), N^ad is the number of hadronic decays, e is the efficiency 
for finding a single lepton from a sample of g^cc Monte Carlo events in which at least one of the 
charmed hadrons decayed semileptonically and passed the selection criteria, and B(c — > Jiiv) is the 
charm hadron semileptonic branching ratio of (9.5 it 0.7) % obtained by taking the average of the 
most recent measurements of OPAL ^ and ARGUS |2l[. With N°^i = 688 ± 42, N^^j = 836 ± 59, 
= (2.72 ± 0.12) % and = (2.83 ± 0.13) % we obtained 



3CC 



0.0309 ± 0.0029, (4) 
: 0.0360 ± 0.0038, (5) 
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where the uncertainties are the statistical contributions only. Figure ^(a) gives the comparison be- 
tween the background subtracted data distribution of the three-jet mass variable, and the JETSET 
prediction for the same distribution (normalised to the data sample). A clear enhancement in the 
mass distribution is visible at low values of the mass ratio, which is well described by the Monte Carlo 
prediction, and when compared to Figure |2|(a), justifies the cut used in section 4.2. Furthermore, 
the agreement between data and Monte Carlo suggests that JETSET describes the production and 
shape of secondary heavy quarks rather well. Figure |^(b) shows the spectrum of the lepton transverse 
momentum for the background subtracted data and for the g^cc Monte Carlo simulation scaled to 
the data size. Here too, a reasonable agreement between the data and Monte Carlo is seen. Leptons 
were accepted if their transverse momentum was smaller than 2.75 GeV/c (where the majority of the 
signal was found). 



5 The D*"^ Analysis 

As D*^ mesons are copiously produced from charm quarks and have a clear signature, they can be 
used to tag charmed hadrons in a gluon jet. Using the three-jet event sample, we chose the gluon jet 
candidate using the jet splitting technique described in section 3.2, and searched for a D*^ meson in this 
jet. The D''+^D°7r+, D° ^K"7r+ decay channel gives a clean signal, since the small mass difference 
between the D*^ and the D*^ limits the phase space available, reducing combinatorial background. 

The D" reconstruction was performed as in |5[ by trying all combinations of oppositely charged 
tracks, assuming one of them to be the kaon, and the other to be the pion. We then added a third 
track, the "slow pion" , demanding its charge to be equal to that of the pion candidate track, to form 
the D*^ candidate. 

The following mass cuts were applied: 

• The reconstructed D*^ mass must lie within 75 MeV of the nominal D*^ mass. 

• The mass difference, AM, between the D*^ and the D" candidates must be in the range 
0.143 GeV/c^ < AM < 0.147 GeV/c^. 

To further reduce the background from random combinations, the following criteria were imposed: 

• The measured rate of energy loss for the kaon candidate track was required to be consistent with 
that expected for a kaon with a probability of more than 0.1. 

• The kaon candidate track momentum satisfied pK > 1.5 GeV/c. 

• The ratio between the measured D*^ energy to the total energy measured in the selected jet 
satisfied Xjet = 1]^^ > 0.2. 

• The helicity angle, 9*, measured between the kaon in the D'^ rest frame and the D'^ direction in 
the laboratory frame satisfied cos 6* < 0.7. The D" decays isotropically in its rest frame creating 
a flat distribution of cos 6*, while background events are peaked at cos 9* = 1. 

Charmed hadrons which pass these selection criteria can also be produced in b hadrons decays. The 



artificial neural network described in Section was used to reject such events. 

Applying all these selection criteria reduced the three-jet data sample from 1.32 million events 
to 308 events. This large reduction in sample size is partially due to the low branching ratio of the 
decay chain as well as the result of the background suppression criteria. The shape of the remaining 
combinatorial background was described by a function of the form (AM — 0.139)"(6-|-cAM-|-(i(AM)^), 
with o, b, c, and d determined from the sideband of the D*^ mass distribution with the signal being 
found in the 2.16 GeV/c2< Mk-^+ < 2.46 GeV/c^ region. The normalisation of the background was 
determined outside the signal region. The estimated background in the signal region was estimated 
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Observed events 


308 


Combinatorial background 
Jet mis-assignment (bb) 
Jet mis-assignment (cc) 
g ^ bb 


171.1±13.1 
19.2±4.3 
35.9±6.0 
4.0±2.0 


Estimated signal 


77.8±15.2 



Table 2: Summary of observed events and estimated background for the D*^ analysis. 



to be 171.1ibl3.1 events. The distribution of AM is shown in Figure ^(a) where a significant signal 
above the fitted background estimation is observed around 145 MeV/c^. 

The jet mis-assignment background was estimated from the Monte Carlo simulation as in 4.3. We 
estimated this background to be 35.9±6.0 and 19.2ib4.3 events from cc and bb events respectively. 

The remaining source of background considered was from secondary b quark pair production 
(g^bb). To estimate its contribution, an average value of g^^^ = (2.69 ± 0.67) x lO^'^ was 
used, along with a selection efficiency of 0.034% (including branching ratios), resulting in a g— >bb 
background estimate of 4.0ib2.0 events. A summary of the sample composition is given in Table |2|. 



In a similar fashion to Section 4.4, we have compared the Monte Carlo simulation with the data. 
Figure |^(b) shows the ratio of the D*^ candidate energy to the jet energy for candidates in the selected 
gluon jet after applying the combinatorial rejection criteria. The Monte Carlo distribution contains a 
g— >cc component scaled to the result obtained for this channel (gcc = 0.04). 

To obtain the rate of g^cc we modified equation 3 to account for the c^D*^, D*^^D'^7r+and 
D'^^K^7r+ branching ratios. The product of the first two was taken to be (15.27ib0.92) % [|l3| and 
B(D°^K-7r+) to be (3.85±0.09) % Q. The number of selected events was N^ei = 77.8 ± 15.2, and 



the selection efficiency found from the Monte Carlo was e = (3.70 it 0.12) %. This gives 

gcc = 0.0408 ± 0.0122, (6) 

where the uncertainty is statistical. 



6 Systematic Uncertainties 

Possible sources of systematic uncertainty and their effect on the background and efficiency are dis- 
cussed below. The systematic uncertainties on gcc from the different sources are summarized in Table |3| 
separately for the electron, muon and D*^ measurements. 

6.1 Lepton Identification 

Muon identification efficiency: The systematic uncertainty from the muon identification efficiency 
was evaluated using a method similar to that in Reference The muon identification efficiency was 
compared in data and Monte Carlo using various control samples, including Z*^— ^/x"*"/!" events, and 
muons reconstructed in jets. Without using dE/dx information, an uncertainty of ±2.0% was found. 

However, as dE/dx information is an important input to muon identification, the effect of mis- 
modelling of the dE/dx was studied. The mean dE/dx for muons in Z*^— events was observed 
to be shifted by approximately 15 % of the dE/dx resolution with respect to the theoretically expected 
value. A similar shift was observed in the Monte Carlo simulation, both for muon pairs and for muons 
identified inclusively in jets, and an uncertainty on the modelling of the dE/dx mean of ±5% was 
assigned. The dE/dx resolution was studied in the data and Monte Carlo using test samples and was 
found to be modelled to better than 5 %. The uncertainties on these sources gave a total uncertainty 
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on the efficiency of muon identffication of ±3.0%. 



Muon mis-identification rate: There was an uncertainty on the number of hadrons that were 
incorrectly identified as muons. This predominantly arises from uncertainties in the simulation of the 
fake probabilities. Three samples were used to test the Monte Carlo modelling of these probabilities: 

• pions from identified Kg^vr+vr" decays. The Monte Carlo predicted that the tracks in this 
sample were more than 99% charged pions with less than 0.1% contamination from prompt 
muons. 

• three-prong r decays. The Monte Carlo predicted that the tracks in this sample were more than 
99% charged pions with less than 0.1% contamination from prompt muons. 

• a sample of tracks passing a set of dE/dx requirements designed to enhance the fraction of 
charged kaons. 

These samples were obtained and used as in Reference [^] to obtain the fake rate uncertainty for 
identified muons from semileptonic decays in Z*^— >bb and Z''— >cc events. This gave a multiplicative 
correction to the Monte Carlo simulation fake rate probability of 1.09 it 0.10. The uncertainty on the 
correction factor strongly depends on the particle composition of the hadron sample. 

This procedure was repeated for the lower momentum spectrum and the particle composition ap- 
plicable for this analysis to give a correction factor of 1.09 it 0.06. As the fraction of pions in this 
analysis is larger than in |18|, the contribution of the relatively large uncertainty associated with the 



charged kaon sample becomes smaller and the overall uncertainty on the fake probability is reduced. 
Hence, the systematic uncertainty on the number of hadrons that were mis-identified as muons is 154 
events. 



Electron identification efficiency: The uncertainty on the simulated electron identification effi- 
ciency was taken from Reference |^^. This study is summarised below. 

The most important variables that were used in the electron identification neural network were the 
specific energy loss dE/dx, its error, and the ratio of the track's energy deposited in the calorimeter 
to the track's momentum. The Monte Carlo simulation and data distributions of these variables were 
compared using samples of identified particles. 

The dE/dx measurements were calibrated in data using samples of inclusive pions at low momenta 
and electrons at 45 GeV/c from Bhabha events. The quality of the calibration was checked with control 
samples, the most important of which were pions from Kg decays and electrons from photons converting 
in the detector. There was a smaller than 5 % difference between the mean dE/dx measured in these 
samples in the data, and the corresponding sample in the Monte Carlo simulation. Similarly, the 
dE/dx resolution in these samples have been studied and the data and Monte Carlo simulation were 
found to agree to within 8%. The uncertainty on the electron identification efficiency from these two 
sources was found by varying both simultaneously, and was ±2.5 % 

A similar study has been performed for the next most significant input variable E/p, which has 
a resolution in the Monte Carlo around 10 % worse than in the data. To correct for this the Monte 
Carlo has been reweighted to match the data, resulting in a variation of the efficiency of ±2.7%. 

No significant contribution to the electron efficiency uncertainty was found from the other input 
variables. The uncertainty arising from them was estimated from the statistical precision of these 
tests, which was less than 1 % of the efficiency. In total, an uncertainty of ±4.0 % was assigned to the 
electron identification efficiency. 

Electron mis- identification rate: The uncertainty on the simulated electron fake rate was eval- 
uated as in [|^]. The study used Kg— >7r'''7r~ decays and three-prong r decays, and gave a fake rate 
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uncertainty of 21 % which corresponds to 17 events. 



Photon conversion tagging efficiency: The uncertainty on the modelhng of the photon conversion 
tagging efficiency was estimated by comparing data and Monte Carlo samples of identified electrons 
with low momentum and low transverse momentum. These samples had a very high electron purity 
and a photon conversion purity of 77%. The uncertainty on the photon conversion tagging efficiency 
was estimated to be 3.5% giving a systematic uncertainty of 27 events. 

Dalitz decays: The uncertainty from this source arises from the modelling of the efficiency and 
from the limited knowledge of the tt^ and ry multiplicities. A systematic uncertainty was assigned 
for possible differences between the multiplicity of neutral pions and rj mesons in the Monte Carlo 
simulation and data. The production rates in the Monte Carlo were varied within the experimental 
uncertainties on the measured multiplicities [Q, giving a 1.2% uncertainty on the electron channel 
result. Combining these sources of uncertainty, we estimated the uncertainty on the number of events 
from Dalitz decays to be 13 events. 

Lepton transverse momentum: To estimate the effect of the lepton transverse momentum criterion 
on the efficiency we compared the fraction of events that pass this criterion in data and in the Monte 
Carlo simulation. We found this effect to be of the order of 1.3%. 



6.2 D*='= Identification 

£)*± reconstruction efficiency: As the efficiency to reconstruct a D*^ was found to be smaller for 
low D*^ momentum, the compatibility of the data and the Monte Carlo D*^ reconstruction efficiency 
was compared by studying the ratio of the D*^ yield in Monte Carlo and data in different momentum 
regions. Since no significant difference between low and high momentum regions was observed, the 
statistical precision of this test was assigned as the systematic uncertainty. The mis-modelling of AM 
resolution was also studied. The observed difference in AM width between the data and Monte Carlo 
was found to be 0.15MeV/c^ resulting in a 2.3% uncertainty on the D*^ efficiency. The uncertainty 
arising from the modelling of dE/dx was estimated by comparing the selection to a modified version, 
where no dE/dx requirements where made. The change in the yield in the Monte Carlo and data 
was consistent. The 3.2% statistical precision of this test was assigned as the uncertainty from this 
source. Combining these effects, we obtained a 4.2 % uncertainty on the D*^ efficiency. 

D*-*- background modelling: The background fitting procedure was repeated changing the region 
of AM used for the normalisation. In addition we also changed the D'' mass sideband to include the 
lower sideband (1.26 GeV/c^< MK-7r+ < 1.56 GeV/c^). Adding the differences due to the above 
tests, and due to the uncertainties on the fitting parameters, we obtained an uncertainty of 3.3 % on 
the number of D*^ mesons. 



6.3 QCD and Fragmentation 

Secondary charm production modelling: As there are no measurements of the momentum spec- 
trum of charmed hadrons from the process of g— >cc to test the modelling, the JETSET parameters 



with which the events were generated were varied according to Reference [22]. The parton shower A 
value was varied in the range 0.13 to 0.31 GeV. The invariant mass cut-off of the parton shower Qo 
was varied between 1.4 and 2.5 GeV. The parameter cjg, the width of the primary hadrons transverse 
momentum distribution was varied in the range 0.37 to 0.43 GeV and b, the parameter of the Lund 
symmetric fragmentation function was varied in the range 0.48 to 0.56GeV-2. Full detector simu- 
lation was not available for all these variations of the model parameters, and therefore, estimates of 
their effect were made by applying appropriate cuts and smearing the event properties. In addition we 
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have compared HERWIG |2| and ARIADNE [||] to JETSET. We have taken the largest difference 
between the three predictions as an extra source of uncertainty. 



Fragmentation modelling: The heavy-quark fragmentation was simulated using the function of 
Peterson et al. and light-quark fragmentation was simulated according to the Lund symmetric 
scheme. The heavy-quark fragmentation model parameters were varied to change the mean scaled 
energy of weakly-decaying bottom and charm hadrons within their experimental range: {xE)h = 
0.702 lb 0.008 and {xe)c = 0.484 it 0.008 respectively In addition, the heavy-quark model was 

changed to the Lund symmetric model, to that suggested by Collins and Spiller and to that of 



Kartvelishvili et al. |26| with parameters tuned according to Reference [27|. The largest difference 
between the Peterson fragmentation and the other models was taken as a systematic uncertainty. 
Jet scheme dependence: The stability of the results has been checked by changing the jet finding 
algorithm to the Durham scheme, with ycut set to ycut = 0.015; this choice of ycut value was made as it 
optimised the significance of the observed signal. With this value, we obtained g^^ = 0.0330 ± 0.0032 
and g[f^ = 0.0401 it 0.0043 for the electron and muon channels respectively (here the uncertainties are 
statistical only). The slightly larger statistical uncertainties with respect to the Jade jet finding results 
(equations 4 and 5) justify the choice of the Jade algorithm for this analysis. Considering the 57% 
overlap of the Jade and Durham samples, the two results are consistent and no additional systematic 
uncertainty was introduced. 

6.4 Heavy Quark Production and Decay 

Semileptonic decay modelling (lepton channels only): Events with a prompt lepton were 
reweighted as a function of the lepton momentum in the rest frame of the decaying heavy hadron 
to simulate different models of semileptonic decay as in [13|. The semileptonic decay model of 



Altarelh et al. |8[ (ACCMM), with parameters tuned to CLEO data |2g] for b decays and to DELCO 
and MARK III data for charm decays, was used for the central values, and was combined with 



the b — > D spectrum measured by CLEO for b ^ c — > £ decays. The model of Isgur et al. [^] 
(ISGW) and their modified model (ISGW**) with the fraction of D** decays determined from CLEO 
data were used to determine the systematic uncertainty due to the b — > ^ spectrum. 



Charm and bottom branching ratios: The dependence on the semileptonic branching ratios, 
b — > £, c — > ^, b ^ c ^ ^, as well as the hadronic branching ratios of b— >D*^, c— >D*^, D^^^D'^vr"'" 
and D'^— >K~7r"*", has been investigated by varying them within their experimental uncertainties [^, 



A possible energy dependence of the B(c^D*^)x B(D*^^D^7r) rate has been taken into account 



by comparing the average of LEP measurements, 0.1527ib0.0092 [13|, with an average of lower energy 



measurements, 0.184 it 0.014 |14|. Typical D*^ energies in this analysis were between these extremes, 
and the uncertainty on this product branching ratio has been inflated to 0.015 to account for these 
possible effects. In contrast, the OPAL and ARGUS [20, ^| values for the semileptonic branching 



ratio c — > ^ are consistent, and no additional uncertainty has been assigned. 

Partial hadronic widths: The partial hadronic widths of bottom and charm with which both the 
efficiency and the jet mis- assignment background were estimated, were taken from a combination of 
LEP and SLD results [O]. The error on the partial widths is a source of systematic uncertainty on 



g — >hh: The uncertainty on the average measured value of g^^|^ gave an uncertainty of 14, 15, and 1.1 
events on the number of background events for the electron, muon, and D*^ channels respectively. 
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Jet mis-assignment background modelling: This background could be mis-modelled if the effi- 
ciency for tagging a heavy quark jet was incorrect, or if the fraction of primary heavy quark jets which 
were identified as the gluon jet candidates was mis-modelled in the Monte Carlo. We tested these 
cases, using enriched Z°— s-bb samples (as described in section 4.4), As data and Monte Carlo were 
consistent in all cases we assigned the statistical precision of these tests as a systematic uncertainty. 
We assigned the difference solely to the signal, using the test purity, and took the larger difference as 
a systematic uncertainty. 



Jet mass criterion: The systematic uncertainty associated with the jet mass criterion was found 
by comparing the fraction of events passing this selection in data and Monte Carlo. The systematic 
uncertainty assigned to this source is 2.3%. 

B hadron decay product multiplicity: The artificial neural network used to reject Z*^^bb events 
was sensitive to the B hadron charged track multiplicity. Consequently Z*^— >bb events were reweighted 
to reproduce the experimental uncertainty on the measured multiplicity |jT^. The uncertainty on the 
number of background events rejected associated with this procedure was estimated to be 5%. 



6.5 Other sources 

Detector modelling: The resolution of the central tracking in the Monte Carlo had an effect on the 
predicted efficiencies and background. The simulated resolutions were varied by ±10% relative to the 
values that optimally describe the data following the studies in [17|. The analysis was repeated and 
the efficiencies and background estimation were recalculated. 



Monte Carlo Statistics: This was the uncertainty due to the finite size of the Monte Carlo samples 
used to determine the efficiencies and background. 



7 Results and Conclusions 

The three measurements of the rate of secondary charm quark production are 

g^_ = 0.0303 ± 0.0028 ± 0.0040, (7) 
g^j. = 0.0353 ± 0.0037 ± 0.0078, (8) 
gD* = 0.0408 ± 0.0122 ± 0.0069, (9) 

where the first uncertainty is statistical and the second systematic. By averaging the leptonic channels 
and taking into account correlations in the systematic uncertainties we obtain 

g^j. = 0.0311 ±0.0022 ±0.0041, (10) 

and combining this with the hadronic channel gives 

gcc = 0.0320 ± 0.0021 ± 0.0038. (11) 

These results can be compared with the previous OPAL measurements [^,|6|. As the values used 
for the rate of the process g— >bb and the branching ratio B(c ^ have changed since the previous 
publication, a comparison should be made after correcting for these effects. The old lepton results |^ 
should be scaled by a factor of 1.09. Thus the scaled results are g^^ = 0.0238 ± 0.0033 ± 0.0045, 
g^j. = 0.0309 ± 0.0063 ± 0.0111 and g^^* = 0.044 ± 0.014 ± 0.015. The analyses described in this 
document use a larger data sample and different selection criteria than before. In particular, the 
different j/cut value and lepton momentum range resulted in a small overlap of the present data sample 
and the previous analysis sample. With these changes, less than one third of the current electron 
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Table 3: Summary of the systematic uncertainties on the measured gcc values. 
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analysis data sample overlaps the previous electron analysis. The increased momentum range of 
muons in this analysis resulted in having only 15 % of this muon sample common to the previous 
analysis. The previous D*^ analysis used only 1.25 million hadronic decays. Therefore, the 
number of common event to the present D*^ analysis is lower than 30 %. Taking into account these 
limited overlaps, the present and previous results are statistically consistent. 

In recent years, the OPAL data calibration has also been refined and the Monte Carlo simulation 
of the detector response had improved significantly. In view of these changes, the analysis presented in 
this paper is more accurate than the previous analysis in both the statistical significance and the proper 
description of the experimental environment. In conclusion, although the two results are statistically 
largely independent, they are dominated by systematic uncertainties, which are better understood for 
the new measurement. This result therefore supersedes the previous OPAL measurement. 

The result presented in this paper, gcc=0.0320ib0.0021ib0.0038, is higher than all the theoretical 
predictions. In particular, the most recent prediction (gcc=2.007% [^), which is higher than most 
other predictions, is 2.7 standard deviations below the result presented in this paper. 
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Figure 1: Jet selection properties: a) number of jets in the events; b) ycut value where a selected jet 
splits into two subjets for all thrcc-jet events; c) jet energy for quark and gluon jets in Monte Carlo 
g— j'ccevents; d) j/cut value at the point where a jet splits into two subjets for quark and gluon jets. 
The Monte Carlo histograms in a) and b) include the default JETSET value of gcc=0.014. 
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Figure 2: Comparison of data and Monte Carlo for the lepton analysis (histograms represent the Monte 
Carlo and data are represented by the points, the dashed histograms show the g^cc spectrum). Monte 
Carlo distributions (solid and dashed histograms) are normalised to the number of hadronic events in 
the data, a) Ratio of the maximum mass of the first two jets with respect to the third jet for events 
that passed the lepton selection and the b-tagging neural network rejection; b) ratio of the multiplicity 
of tracks and unassociated clusters in the first two jets with respect to the third jet for the events 
in a); c) lepton candidate spectrum for events passing all the selection criteria; d) lepton candidate 
transverse momentum spectrum for events passing all but the transverse momentum criteria with the 
g^cc spectrum scaled to the signal area. The arrows in a), b) and d) show the cut value used in this 
analysis. Events below these values were accepted. 
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Figure 3: Distributions of a) Max(Mi, M2)/M3, b) lepton transverse momentum for background sub- 
tracted data (points) and for g— >cc Monte Carlo scaled to the average lepton analysis gcc result (solid 
line): All the selection criteria described in the text were applied with the exception of the variable 
in each plot where the cut is shown by the arrow. The error bars shown represent the statistical and 
systematical uncertainties combined. 
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Figure 4: (a) The distribution of the difference of the invariant mass between the D* candidate and 
the D'^ candidate, AAf. The sohd curve represents the estimated background shape. An estimated 
137±13 D*^ candidates are above the background fit around 0.145 GeV/c^. (b) Ratio of the D*^ can- 
didate energy to the jet energy for candidates in the selected gluon jet after applying the combinatorial 
rejection criteria. The solid line shows the Monte Carlo spectrum with the g^cc component scaled 
to measured value in the D*^ analysis (gcc=0.04), and the dashed histogram shows the spectrum for 
true D*^ mesons from g— >cc scaled to the measurement area. 
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